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Optically-induced nuclear spin polarization in a fluorine-doped ZnSe epilayer is studied by time- 
resolved Kerr rotation using resonant excitation of donor-bound excitons. Excitation with helicity- 
modulated laser pulses results in a transverse nuclear spin polarization, which is detected as a change 
of the Larmor precession frequency of the donor-bound electron spins. The frequency shift in depen¬ 
dence on the transverse magnetic field exhibits a pronounced dispersion-like shape with resonances 
at the fields of nuclear magnetic resonance of the constituent zinc and selenium isotopes. It is 
studied as a function of external parameters, particularly of constant and radio frequency external 
magnetic fields. The width of the resonance and its shape indicate a strong spatial inhomogeneity of 
the nuclear spin polarization in the vicinity of a fluorine donor. A mechanism of optically-induced 
nuclear spin polarization is suggested based on the concept of resonant nuclear spin cooling driven 
by the inhomogeneous Knight field of the donor-bound electron. 

PACS numbers: 76.60.-k, 76.70.Hb, 78.47.D-, 78.66.Hf 


I. INTRODUCTION 

The hyperfine interaction between electron and nuclear 
spins in semiconductor structures has been of particular 
interest over many years [I|, 0. Lately it was intensively 
driven by the intention to extend the electron spin co¬ 
herence time and the idea to employ the nuclear spins 
as quantum bits (qubits) Q. Using optical excitation 
with circularly polarized light angular momentum can 
be transferred via the electron system to the nuclei. The 
polarized nuclei, in turn, act back on the electrons as an 
effective magnetic field (the Overhauser field), causing a 
splitting of the electron spin states [l|. Without an ex¬ 
ternal magnetic field, the hyperfine interaction with the 
nuclear spins is the main source of dephasing for local¬ 
ized electron spins in which fluctuations in the nuclear 
spin polarization induce fluctuations in the electron spin 
splittings d^. 

The nuclei-induced electron spin splitting can reach 
values of 100 ^eV for high nuclear spin polarization. It 
can be detected spectrally as a splitting of the emission 
line, if the linewidth of the photoluminescence (PL) is 
sufficiently narrow, as it is the case for single emitters 
like quantum dots (QDs) or impurity centers. A rela¬ 
tively high nuclear spin polarization can be created using 
the method of optical pumping in longitudinal magnetic 
field (parallel to the optical axis and orthogonal to the 
sample surface) In the case of inhomogeneously 

broadened optical transitions, as it is the case in emit¬ 
ter ensembles, optical methods do not provide sufficient 
spectral resolution so that alternative approaches should 
be used. The nuclear spin polarization can be indirectly 
measured by its influence on the electron spin polariza¬ 
tion in a transverse magnetic field, the Hanle effect [HI- 


In this case the Overhauser field modifies the measured 
degree of the circular PL polarization as a function of 
magnetic field applied in the direction transverse to the 
optical excitation [3, . If one additionally applies a 

polarization modulation of the excitation at a frequency 
corresponding to the nuclear magnetic resonance of one 
of the constituent isotopes, a nuclear spin polarization 
in the direction transverse to the optical excitation can 
be achieved [iMl- However, such measurements are 
limited to the field ranges within the width of the Hanle 
curve, which is determined by the spin dephasing time 
of the carriers, their g factor values and the nuclear spin 
polarization. 

In the recent publication of Zhukov et al. the au¬ 
thors extended the magnetic field range by implementing 
a novel experimental technique to measure the optically- 
induced nuclear magnetic resonances by assessing the 
electron spin coherence in the resonant spin amplification 
(RSA) regime. For that purpose time-resolved pump- 
probe Kerr rotation (TRKR) was used [U, [2^ . Informa¬ 
tion about the Overhauser field that changes the electron 
spin splitting was thereby transferred from the spectral 
domain into the temporal domain. In the present paper 
we extend the investigation of dynamic nuclear spin po¬ 
larization (DNP) by applying the developed experimen¬ 
tal approach to ZnSe doped with fluorine donors. We 
observe changes of the electron Larmor precession in a 
broad range of transverse magnetic fields, which depends 
on the pump power. The induced shifts of the Larmor 
precession change sign at the fields of the nuclear mag¬ 
netic resonances (NMRs). A theoretical approach to an¬ 
alyze the mechanism of dynamical nuclear spin polariza¬ 
tion for our experimental conditions is developed. 

Fluorine donors in ZnSe (further ZnSe:F) and the cor- 
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responding donor-bound electron spins are currently of 
substantial interest because of several recent key achieve¬ 
ments towards solid-state quantum devices: Sources of 
indistinguishable single photons [2^ as well as entangled 
photon-pa.irs and optically controllable electron-spin 
qubits |25l - l^ have been demonstrated so far. The spin 
coherence of the donor-bound electron is generally lim¬ 
ited by the non-zero nuclear spin background in the host 
crystal Q. However, in ZnSe isotopic purification can 
be applied to deplete the remaining low amount of non¬ 
zero nuclear spins. Apart from that, fluorine has a nat¬ 
ural 100% abundance of spin 1/2 nuclei, which might be 
considered as a nuclear spin qubit coupled to a single 
electron spin qubit via the hyperfine interaction. These 
aspects make the ZnSe:F system particularly attractive 
for the investigation of electron and nuclear spin related 
features. 


II. EXPERIMENTAL DETAILS 

The sample under study is a homogeneously fluorine- 
doped, 70-nm-thick ZnSe:F epilayer grown by molecular- 
beam epitaxy on a (OOl)-oriented GaAs substrate. The 
epilayer is separated from the substrate by a 20-nm-thick 
Zno.ssMgg j^gSe barrier layer to prevent carrier diffusion 
into the substrate. The barrier in turn is grown on top of 
a thin ZnSe buffer layer to reduce the defect density at 
the III-V/II-VI heterointerface. The concentration of flu¬ 
orine donors in the ZnSe:F epilayer is about 10^®cm“^, 
so that the distance between the neighboring donors is 
larger than the Bohr radius of the donor-bound elec¬ 
trons The sample was placed in an optical cryostat 
with a superconducting split-coil magnet. The sample 
temperature was fixed at T = 1.8 K in all measurements. 

Figure [TJa) shows the normalized PL spectrum mea¬ 
sured for continuous-wave (cw) excitation with a photon 
energy of 3.05 eV. The PL is detected with a Si-based 
charge-coupled-device camera attached to a 0.5 m spec¬ 
trometer. In the characteristic peak pattern each feature 
can be assigned to different exciton complexes [2^. The 
labels in the figure mark the following optical transitions: 
FX refers to the free exciton and D°X to the donor-bound 
exciton containing heavy holes (HH) or light holes (LH), 
respectively. 

To obtain insight into the electron spin dynamics we 
use the TRKR technique. The electron spin coherence is 
generated by circularly-polarized pump pulses of 1.5 ps 
duration (spectral width of about 1 meV) generated by 
a mode-locked Ti:Sa laser operating at a repetition fre¬ 
quency of 75.75MHz (repetition period Tr = 13.2ns). 
The laser frequency is doubled by a BBO (beta barium 
borate) crystal to convert the Ti:Sa range of photon en¬ 
ergies from about 1.25 — 1.7eV to 2.5 — 3.4eV. After 
excitation of the sample along the growth axis z with 
the circularly-polarized pump pulses, the reflection of the 
linearly-polarized probe pulses is analyzed with respect 
to the angle of polarization rotation as a function of delay 
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FIG. 1. (a) Normalized PL spectrum measured at R = OT 
and T — 1.8 K. (b) Kerr rotation signal measured for resonant 
D'^X-HH excitation at 2.800 eV. B — 0.5 T. The arrow marks 
the negative time delay of the probe with respect to the pump 
at which the RSA signal is measured. 


between the pump and probe pulses. The pump helicity 
is modulated between cr+ and cr~ circular polarization 
by an electro-optical modulator (EOM) using modula¬ 
tion frequencies /m in the range 20 — 500 kHz. The probe 
beam is kept unmodulated. The Kerr rotation (KR) an¬ 
gle is measured with a balanced detector connected to a 
lock-in amplifier. 


Figure [TJb) shows a TRKR signal measured at a mag¬ 
netic field of B = 0.5 T applied perpendicular to the 
growth axis z (Voigt geometry). Both pump and probe 
have the same photon energy (degenerate pump-probe 
scheme) and are resonant with the D°X-HH transition 
at 2.800eV. The pump power is P = I.lmW and the 
probe power 0.5 mW, both having a spot size of about 
300 /rm on the sample. The oscillations in the TRKR 
signal correspond to the Larmor spin precession of the 
donor-bound electron with a g factor of ge = 1-13 ± 0.02. 
The oscillating signal at negative pump-probe delays in¬ 
dicates that the spin coherence is maintained up to the 
arrival of the next pump pulse, i.e. the electron spin de¬ 
phasing time T 2 * is close to the laser repetition period Tr. 
Due to the long spin dephasing time the action of the suc¬ 
ceeding pump pulses on the remaining spin polarization 
depends on the electron g factor and the magnetic field 
strength. In this case, the resonant spin amplification 
regime can be used for studying spin coherence. 

To that end we hx the probe at a slightly negative time 
delay of At = —20ps with respect to the pump pulse 
arrival moment (see arrow in Fig. [Ijb)) and detect the 
KR signal while scanning the transverse magnetic field. 
Further, we conducted experiments where an additional 
oscillating magnetic field (radio frequency or RF-field) 
Brf = (0, 0,i?RF,z) is applied along the z axis using a 
small coil placed close to the sample surface. 
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FIG. 2. (Color online) (a) Normalized RSA signals measured 
at different pump powers and fm = 185 kHz. The NMR of 
^^Se at B = 22.76 mT is marked by the arrow. The left inset 
shows a close up of the RSA peaks around 29 mT. The posi¬ 
tion of the peak maximum is plotted as function of the pump 
power in the right inset. T — 1.8 K. (b) Black solid circles 
show how the difference (shift) between the peak positions at 
8mW and at 0.1 mW from panel (a) changes in dependence 
on the magnetic field. Open circles demonstrate the induced 
shift measured directly by time-resolved KR (see text). Lines 
are guides to the eye. 


III. EXPERIMENTAL RESULTS 


A. Observation of the nuclear spin polarization 

Figure ura) shows a set of RSA curves of normal¬ 
ized amplitude for an EOM modulation frequency of 
/ni = 185 kHz, measured for different pump powers. The 
RSA curves have the typical shape of periodic peaks as 
function of the magnetic field strength. The dependen¬ 
cies of the RSA peak distance and width on magnetic 
field are determined by the electron g factor the g^ 
spread Ag, the spin dephasing time T, and the repeti¬ 
tion period of the laser pulses Tr [^,[3. Further optical 
properties of this sample and information on the electron 
spin dephasing and relaxation mechanisms can be found 
in Refs. [Hlail. 

Let us turn now to the experimental results on nuclear 
effects induced and detected via the spin dynamics of the 
resident donor-bound electrons. As mentioned above, 
the sample is excited by the helicity modulated pump. 
Under these conditions nuclear spin polarization in the 
system is usually suppressed, as the nuclei cannot fol¬ 
low the oscillating electron spin polarization fast enough 


due to the long nuclear spin relaxation time. However, 
nuclear spin polarization can be observed for finite mag¬ 
netic fields in close vicinity to the NMR frequencies (cor¬ 
responding to the helicity modulation frequency) of the 
constituent isotopes. Here the nuclear spin system be¬ 
comes polarized through its cooling in the magnetic field 
which results from the electron spin polarization (Knight 
field) [^ . 1^ . The induced nuclear spin polarization ap¬ 
pears in RSA curves as an additional amplitude mod¬ 
ulation at the NMR fields, see in Fig. [2Ja) the narrow 
resonance for the ^^Se isotope at 22.76 mT measured for 
/m = 185 kHz at pump powers above 0.5 mW. 

If we consider the behavior of the RSA peaks at dif¬ 
ferent excitation powers we observe the following behav¬ 
ior: Fig.[2Ua) displays several RSA curves taken at pump 
powers varied from P = 0.1 mW (bottom spectrum) up 
to 8mW (top spectrum). The probe power is kept in 
the range of 0.5 mW for all measurements. The signal 
amplitudes are normalized and the curves are displaced 
vertically relative to each other, to simplify comparison 
of the peak positions. The increase of the pump power 
leads to a shift of the RSA peaks, which reflects a change 
of the electron precession frequency. The left inset in 
Fig. Ha) provides a closeup of a single peak to highlight 
this shift. The peak positions for each curve are given 
in the right inset, which evidences a saturation effect at 
high pump powers. 

One can see from Fig. Ha), that the direction and the 
magnitude of the shift depend on the RSA peak position 
relative to the NMR field and also on the pump power. 
The peaks located at lower fields relative to the resonance 
exhibit a shift towards lower magnetic fields for increas¬ 
ing pump power, while the peaks at higher fields than 
the resonance are shifted towards higher magnetic field. 
The shift itself can be explained by an additional induced 
magnetic field acting in addition to the external mag¬ 
netic field on the electron spins. The shift of RSA peaks 
to lower magnetic fields indicates an additional magnetic 
field pointing in the same direction as the external field. 
Vice versa, the shift to higher fields indicates opposite 
orientations of the induced and external fields. 

To demonstrate how the peaks shift across the whole 
range of external magnetic fields, we plot the differ¬ 
ence of the peak position taken at maximal (saturated, 
P = 8mW) and minimal used pump powers (0.1 mW) 
as a function of magnetic field. The black solid circles in 
Fig.Hb) show this dependence. The circles are placed at 
the fields of the RSA peaks measured at minimal pump 
power. The induced RSA peak shift (vertical axis in 
Fig.Hb)) directly corresponds to the opposite strength 
of the induced nuclear field. 

Figure H^-) shows the pump induced shift over a 
broader range of magnetic field, measured at a modu¬ 
lation frequency of /m = 250 kHz. Here the NMR of the 
^^Se and ®^Zn isotopes are observed at \B\ = 30.75 mT 
and 93.75 inT, respectively. FigureHb) shows in addition 
the RSA shifts in the vicinity of the ^^Se NMR measured 
for various /m up to 1 MHz. The amplitude and the shape 
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FIG. 3. (Color online) (a) Induced RSA shift measured at 
/m = 250 kHz in a broader range of magnetic fields compared 
to Fig. [2l Therefore the NMR of the isotope is also 
observed, (b) Shift of the RSA peaks measured at different 
modulation frequencies /m. 


of the shifts remain unchanged in this field range. 

The dependencies in Figs. O^b) and [3] exhibit several 
peculiarities: 

(i) They show a characteristic dispersive shape around 
the NMRs of the constituent isotopes. In Fig. EJa) 
one sees these resonance-like dispersive features for 
both the selenium and zinc isotopes. In the follow¬ 
ing we will call them resonances for simplicity. 

(ii) The resonances occur at magnetic fields, which are 
much larger than the width of the Hanle curve. 
Here the Hanle-curve width is given by the width 
of the RSA peak around zero magnetic field and is 
about 2mT (ful l width half maximum - FWHM), 
see also Ref. (^ . 

(hi) The resonances are extended over a quite broad 
magnetic field range as result of their slowly decay¬ 
ing dispersive tails. 

For a deeper understanding, the open circles in 
Fig. [2jb) show the results of additional measurements, 
where the induced nuclear fields are not only given at the 
RSA peak positions, but are also determined in between 
the peaks to obtain a higher resolution on magnetic field. 
To provide these data we determine the electron preces¬ 
sion frequency from a measurement of the KR signal as 
a function of the time delay between pump and probe 
pulses (see Fig. [Ifb)). We perform the measurements 
for transverse magnetic fields varied from 10 mT up to 
40 mT using a 1 mT incremental step. Each data point 
represents the difference of the electron Larmor frequen¬ 
cies measured at 8inW and O.ImW pump power. Using 
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FIG. 4. (Golor online) (a) Normalized RSA spectra for differ¬ 
ent circular polarization degrees of the pump. The spectra are 
shifted vertically relative to each other for clarity. T = 1.8 K, 
P = 8 mW. (b) Position of the RSA peak around 28 mT held 
strength in dependence on the degree of circular polariza¬ 
tion. Red line is a parabolic ht to the data demonstrating 
the quadratic dependence of the nuclear spin polarization on 
the circular polarization, (c) Relative peak shifts for different 
circular polarization degrees in dependence on the magnetic 
held. The data are normalized to the shift of the peak at 
24 mT. 


the electron g factor = 1.13, one can convert this Lar¬ 
mor frequency difference to a magnetic field, which cor¬ 
responds to the RSA peak shift, so that we can present 
it on the same scale as the RSA measurements. 

One sees that the experimental data given by the open 
circles in Fig. Elb) follow closely the solid circles. Both 
curves show a slight vertical shift to positive values. The 
FWHM of the broad dispersive resonance is on the or¬ 
der of 40 mT. The high resolution of the open circles 
data is particularly advantageous for the transition re¬ 
gion slightly above 20 mT where the sign reversal of the 
RSA peak shift takes place. This sign reversal occurs in 
a field range of 1 — 2 mT only. Note that the long tails 
and the sharp transition cannot be explained by a sim¬ 
ple dispersive curve with a large homogeneous linewidth, 
which will be addressed in more detail by our model con¬ 
siderations in Sec. IIV Al Additionally, one can see that 
the open circles data exhibit oscillations with a period 
equal to the period of the RSA peaks. These oscillations 
could not be resolved in the data set given by the closed 
circles, as there the measurements were performed only 
at the RSA peak positions. 


B. Electron spin polarization 

Before we turn our attention to the investigation of 
the peculiarities listed in the previous section, we check 
how the induced nuclear spin polarization depends on 
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the electron spin polarization. Figure HJa) shows the 
measured RSA signals for several degrees of pump cir¬ 
cular polarization at /m = 185 kHz and P = 8 mW. The 
signals are normalized in amplitude to simplify compari¬ 
son of the peak shift. One clearly sees for the RSA peak 
close to the NMR, that the pump-induced shift decreases 
for lower circular polarization degrees. This confirms 
that the nuclear spin polarization is induced by spin- 
oriented electrons. Figure SJ^b) shows an example of the 
peak shift dependence on the polarization degree around 
28 mT. The red line represents a parabolic fit to the data, 
demonstrating the quadratic dependence of the nuclear 
spin polarization on the degree of circular polarization. 

Additionally, Fig. |4l^c) shows how the peak-shift de¬ 
pends on the magnetic field for three different degrees 
of polarization. Obviously the width of the resonance 
is proportional to the induced nuclear spin polarization 
and decreases for lower polarization degrees of the elec¬ 
tron spin. Here the shifts are given relative to the peak 
position for 15% polarization degree. The shift ampli¬ 
tudes are then normalized to the shift of the peak at 
about 24 mT. This allows us to clearly resolved the ac¬ 
celerated decay with magnetic field for lower degrees of 
the nuclear spin polarization. 

For a complete understanding of the nature of the nu¬ 
clear spin polarization produced by the polarized elec¬ 
tron spins, information about the average electron spin 
polarization S is important (^ . To decide which spin 
component is responsible for the nuclear spin polariza¬ 
tion we should provide a tomographic measurement for 
all three of them: Sx, Sy and Sz- 


1. Sz-component: Time-resolved Kerr rotation 

In our experiment the average Sz component, which 
is created along the optical excitation axis can be evalu¬ 
ated directly from the experimental data by integrating 
the KR signal over the whole time period between the 
pump pulses. This averaging is expected to result in a 
finite value of Sz in magnetic fields for which the Larmor 
precession period, Tl, is longer than or comparable with 
the spin dephasing time T^. this relation is not fulfilled 
in strong magnetic fields, where Tl <C , Sz ~ 0, but in 
the field range studied here it is perfectly valid. 

Figure [SJa) demonstrates an example of the measured 
Kerr rotation signal in the pump-probe delay range from 
0 to 13.2ns at B = 25.5mT (black curve). To evaluate 
Sz, which is directly proportional to the KR amplitude, 
we fit the data with an exponentially decaying cosine 
function and integrate the fit curve over the whole period 
of Tr = 13.2 ns. The magnetic field dependence of Sz is 
given in Fig. [5](b). It has a finite value for B < 15 mT 
and approaches zero for higher fields. It is instructive to 
compare the Sz{B) dependence with the results of the 
RSA shifts shown again in Fig. [51(c). From this compar¬ 
ison one can conclude, that the high nuclear spin polar¬ 
ization, corresponding to large shifts of RSA peaks, is 
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FIG. 5. (Color online) (a) Example of a time resolved pump- 
probe Kerr signal measured at B = 25.5mT (black curve). 
Red curve is a fit to the data with an exponentially decaying 
cosine function recorded over the full range of delays of Tr = 
13.2 ns between two pump pulses, (b) Average electron spin 
polarization along the z axis resulting from the fitted curves 
as a function of magnetic field, (c) Reproduction of the data 
from Fig.l^b) for simplified comparison with the data in panel 
(b). 

induced in magnetic fields, where Sz is already zero. 


2. Sx: Knight field influence 

The presence of an average electron spin polariza¬ 
tion has an effect on the nuclei through providing an 
effective mag netic field, the so called Knight field, see 
Sec. IIV Al [y. To evaluate its value along the external 
magnetic field, which is proportional to Sx spin compo¬ 
nent, one can evaluate the effect of the Knight field from 
the NMR frequency dependence on the magnetic field. 
For this purpose, we have scanned the RF-field frequency 
for C/rf = 0.05 V around the ^^Se NMR for different 
magnetic fields. Any Knight field component produced 
by an average Sx component would induce an additional 
magnetic field along the x-axis and thereby lead to an 
offset of the linear dependence of NMR resonance of the 
^^Se isotope on the magnetic field. Figure IHKa) shows 
the RSA curve measured for /m = 50 kHz at T = 8 mW, 
where the ^’’Se NMR is seen at about Rnmr = 6 mT. Fig¬ 
ure IHKb) gives an example of the RF-frequency scan at a 
fixed magnetic field of R = 5.1 mT. As shown by the red 
fit using a Lorentz curve, the central NMR frequency is 
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FIG. 6. (Color online) (a) RSA curve measured for/m = 
50 kHz at P = 8 mW. (b) RF frequency scan at fixed 
B = 5.1 mT. fc gives the central frequency defined by the 
Lorentzian fit (red curve) to the data, (c) Dependence of the 
fc on the magnetic field. Red line is a linear fit to the data. 

located at fc = 39.37 ± 0.02 kHz, showing also the preci¬ 
sion of the measurement of the NMR position: 0.02 kHz 
accuracy corresponds to about 3 field strength. Fi¬ 
nally, Fig. Uc) represents a collection of all measured fc 
at different magnetic fields. This data set demonstrates a 
linear dependence of the resonance frequency with mag¬ 
netic field where the corresponding fit leading to an offset 
of —0.4±0.2kHz, which corresponds to —51±26^T. This 
allows us to conclude, that the Knight field produced by 
a possible Sx spin component, if present at all, should 
not be the reason for the measured nuclear spin polariza¬ 
tion. As will become clear from the next sections, such 
an offset (or a Knight field) is negligible compared to the 
Knight field produced along the y direction. Namely, at a 
distance of one localization radius a; from the donor cen¬ 
ter, the Knight field reaches a strength of 3.5 mT along 
y direction, see Sec. IIVBI 

3. Sy: RF-field versus Knight field 

To test the presence of a Sy spin component we applied 
a RF-field with the same frequency as the used helicity 
modulation frequency /m. The relative phase between 
this RF-field and the helicity modulation could be con¬ 
trolled, as well as the RF amplitude. 

In the frame system rotating with frequency /m, one 
can interpret the RF-field as an additional, constant-in¬ 
time magnetic field acting on the nuclei (see Sec. IIV Al 
on the rotating frame system). This allow us to find the 
orientation of the average electron spin component and 
compensate its action on the nuclei by counteracting with 
the applied RF-field, see Fig.[7Dc). 

Figure 0 demonstrates the results of such measure¬ 
ments at different magnetic fields for /m = 185 kHz. The 
magnetic field values are comparable to those applied in 


the RSA studies Fig.[2Da). Here we plot the signal deter¬ 
mined as the difference between the KR amplitudes with 
and without the RF field, applied at a fixed magnetic 
field for a high pump power of P = 8 mW. Figure [TKa) 
shows the influence of the RF field on the KR amplitude 
at the field position oi B = 26 mT on the right hand 
side of a RSA peak, i.e. in a field range where the KR 
signal has a decreasing slope with increasing field. The 
color coding is such that blue color means reduction of 
the Kerr rotation signal under the influence of the ap¬ 
plied RF-field, while red color represents an increase of 
the signal. The panels right next to the contour plot 
show cuts through the contour plot, namely the RF am¬ 
plitude dependence at fixed phase of —90° (top) and the 
phase dependence at the saturated RF amplitude level of 
Prf = 5V, which corresponds to an effective magnetic 
field of 170 ^T (bottom) [3^. 

As one can see from the RF dependence, the Kerr am¬ 
plitude is reduced for higher RF voltages, i.e. the RSA 
peak shifts to lower magnetic fields so that the nuclear 
spin polarization is reduced. This is most efficient for 
the phase close to —90°, implying that the Knight field 
component should initially be oriented along the y axis. 
In the transverse magnetic field it becomes then rotated 
to the direction parallel to the z axis, thereby accumulat¬ 
ing the phase of —90°, which is best compensated when 
the RF field, applied along the — z axis, is acting di¬ 
rectly against it, see Fig.[7Kc). As an additional support 
for this interpretation, we observe an increase of the KR 
amplitude at -1-90°, where the nuclear spin polarization 
becomes amplified by the RF field. Here the Knight field 
is acting in the same direction as the RF, see Fig. [7Kc), 
increasing therefore the overall nuclear spin projection 
/N,a; along the i-axis. 

A similar behavior is observed at a higher magnetic 
field strength, B = 28.7 mT, as shown in Fig.jTJb). Here 
we are on the left increasing slope of the next RSA peak, 
so that the overall amplitude changes are inverted. The 
RSA peak shifts to lower magnetic fields for reduced nu¬ 
clear field and the KR amplitude increases. The relative 
phase behavior of the RF however stays the same: we re¬ 
duce the nuclear spin polarization at —90° and increase 
it at -1-90°, which shows that the Knight field has a con¬ 
stant phase shift relative to the RF field and is oriented 
along the ?/ axis. 

Several more RSA peaks at higher magnetic fields were 
tested demonstrating similar behaviors (not shown here). 
Measurements at different magnetic fields demonstrate 
that the electron spin orientation stays the same. Addi¬ 
tionally, the effect of the RF field at the phase of —90° 
leads to similar changes of the KR amplitude in the range 
of fields 20 — 40 mT which has an value of about 1.9 ± 0.3 
(arb. units of KR amplitude). This means, that there 
is a finite average electron spin component along the Sy 
direction, which stays constant when varying the trans¬ 
verse magnetic field strength. 
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FIG. 7. (Color online) (a) Contour plot for the phase between the RF field applied along the z axis and the helicity modulation 
for different RF amplitudes at a magnetic field strength of 26 mT. On the righthand side of the plot cuts through the contour 
at the phase of —90° (top) and at the RF voltage of 5 V (bottom) are shown, (b) Same as in (a) but for a different magnetic 
field position, namely 28.7 mT. /m = 185 kHz. (c) Schemes of the rotating frame systems (RFS) for the RF field acting along 
the Knight field +90° and antiparallel to it —90°. Both schemes are given for B > Hnmr. 


IV. THEORETICAL MODEL AND DISCUSSION 
A. Classical treatment 


Let us consider excitation with circular polarization 
alternated between left and right at the modulation fre¬ 
quency /m. An external magnetic field is applied perpen¬ 
dicular to the pump light k vector, kpump, which is shine 
in parallel to the z axis. Figure [5Ka) shows this configu¬ 
ration in the laboratory frame system (LFS). It is known 
that under helicity modulated excitation with modula¬ 
tion frequency 27r/m an optical polarization 

of the nuclear spins occurs only for B close to the reso¬ 
nance field Bnmr = STr/m/y. Here is the nuclear 

spin relaxation time and 7 is the gyromagnetic ratio of 
the nuclear isotope Measurements of will 

be presented elsewhere. For the studied ZnSe:F epilayer 
these times fall in the range of tens of milliseconds. The 
literature values of 7 for the ^^Se (yse = 51.08 x 10® [T 
s]“^) and ®^Zn isotopes {'^zn = 16.77 x 10® [T s]“^) lead 
to Hnmr [niT]= 0.123/m [kHz] and 0.375/m [kHz], respec¬ 
tively |34l ] . These values are in very good agreement with 
the values observed for the NMR resonances. 

In the following we consider only the ^^Se isotope hav¬ 
ing nuclear spin 7 = 1/2 with an abundance x = 0.0758. 
The hyperfine constant Ase = 33.6 /reV, which is taken 
for a primitive cell with two nuclei |35l . 1361] . 

The dynamic nuclear spin polarization is caused by nu¬ 
clear spin flips in the presence of the electron Knight field 
Bk = bgS which precesses synchronously with the elec¬ 
tron spin S and provides a temporally constant energy 
flow into the nuclear spin system [ 3 ^ . Here 


AseVp 

XSeflTTaf ’ 


( 1 ) 


volume with a two-atom basis and oq = 0.566 nm is the 
lattice constant of ZnSe [s^; ai gives the localization 
radius of an electron at the donor. 

The average nuclear spin polarization In has a compo¬ 
nent along the external field B, /n.s:, and, therefore, af¬ 
fects the electron spin precession frequency, see Fig.[51[b). 
If the external magnetic field is close to the Hnmr, the 
projection of the nuclear field on the external field B is 
given by: 

^ _ (S ■ Bk)(B - .Bnmr) 

where = 0.006 mT is the root mean square local held 
due to the nuclear dipole-dipole interactions 

This can be interpreted as nuclear spin cooling in the 
rotating frame system (RFS) [13, see Fig. [8][b). The 
Knight held, oscillating with frequency 27r/m, can be de¬ 
scribed by a superposition of two helds, each rotating 
around the external held opposite to each other. Close 
to the NMR frequency the component, which rotates in 
the same direction as the nuclear spin is the important 
one, while the other one can be neglected. In ZnSe all nu¬ 
clei have 7 > 0 , therefore this component is rotating 

counterclockwise if one looks in the direction along the 
external held B. In the RFS the electron spin is constant 
and the nuclei see the total magnetic held BTotai, com¬ 
posed of the constant Knight held Bk and the effective 
external held B — Bnmr. This situation is analogous to 
the cooling in the laboratory frame system (LFS) in the 
stationary case. One can see from Eq. that upon 
fulhlling the resonance condition B = T^nmr the nu¬ 
clear spin polarization along the external magnetic held 
7n,x — 0- 

T he p olarized nuclei create an Overhauser held (see 
Ref. [^) which acts on the electron spins: 


characterizes the maximal Knight held amplitude at the 
center of the donor 37|. vq = a^/d is the primitive cell 




AseXlN,x 


(3) 


MB5e 




























FIG. 8. (Color online) (a) Spin orientation in the laboratory 
frame system. The red and blue arrows symbolize the elec¬ 
tron spins generated with cr'*’ and cr“ photons, respectively, 
alternating in helicity with fm- (b) Scheme of magnetic fields, 
acting on the nuclei in the rotating frame. The frame is ro¬ 
tating about the x axis with frequency 2%/^ in the direction 
of the nuclear spin precession. Bk and S are the average 
electron Knight held and spin, respectively, which are hxed 
in the rotating frame and are rotating with 2'Kfni in the lab¬ 
oratory frame, /n.s; is the projection of the induced nuclear 
spin polarization onto the direction of the external held B. 

Here /tb is the Bohr magneton. From Eq. m one can 
see that the following has a dispersive 

shape. Its sign is defined by the detuning H—Hnmr- This 

describes exactly the behavior of the signal described in 
Sec. mill item (i). 

By adding to or subtracting from the external magnetic 
field, the field B^^x alters the electron spin precession 
frequency. In turn, this leads to a change of the elec¬ 
tron spin polarization value at a specific external mag¬ 
netic field. Usually, this effect occurs if the resonant held 
^NMR does not exceed the half width of the electron spin 
depolarization curve i?i /2 (Hanle curve) [13 111- Other¬ 
wise, the electrons become depolarized due to their spin 
precession about the transverse held B -|- Bn,x, so that 
the value (S • Bk) and the resulting Hn x are small. Nev¬ 
ertheless, it has been shown in Ref. [1^ that such reso¬ 
nances do occur also at the Hnmr helds, which are much 
stronger than i?i /2 and have a width of about one mT, 
see item (ii). In what follows, we concentrate on the na¬ 
ture of the broad resonance (item (hi)). 

As one can see, Eqs. m and (|^ support the measure¬ 
ments presented in Figs.HKa) and|3Kb). The nuclear spin 
polarization, which causes the RSA shift via the Over- 
hauser field indeed follows the dependence. The 
5^ dependence comes from the product (S • Bk), where 
Bk = beS. The accelerated decay of the nuclear spin 
polarization with increased magnetic held is also well 
described by the dependence in the denominator of 
Eq. (HI, see Fig. UKc). 

As mentioned before, the sign of the induced shifts in 
Fig. [2{b) is opposite to the induced nuclear Helds, 
Taking into account that in ZnSe the electron g factor 
ge > 0, so that < 0 and that the hyperhne constant 
Ase > 0, Eq. (131) reproduces this sign dependence. It is 
negative for B > Hnmr and positive vice versa. 

We now try to simulate on the basis of Eq. the 


FIG. 9. (Golor online) Peak shift and corresponding Bn, a: 
calculated for the ^^Se isotope after Eq. ([5]) (black line) with 
be = —370mT, S = 0.07 and ai = 3.4 nm. For comparison we 
show two homogeneously broadened curves calculated after 
Eq. (O using fee = — lOmT (red line) and fee = — lOOmT 
(blue line). Open circles represent the experimental data as 
in Fig. HDb), but shifted slightly downwards to compensate 
for the offset. 


data shown by the open circles in Fig. [2Kb) and repro¬ 
duced in Fig.[9l The Htting parameters in this procedure 
are feg and S. S is Hxed at 0.07, which will be JustiHed 
below. Figure O demonstrates examples for a homoge¬ 
neous nuclear spin polarization using fee = — lOmT and 
— lOOmT, given by the red and blue curves, respectively. 
As can be seen, the fee = — lOmT curve describes only 
the fast changeover of the sign close to the resonance 
very well, but completely fails to Ht the data away from 
the resonance. On the other hand, the fee = —lOOmT 
curve shows the right tendency compared to the data 
only in the tails, far from the resonance condition. The 
amplitudes here, however, deviate considerably from the 
measured data points. 

This simulation leads to an important conclusion: the 
Knight Held is not constant within the localization area 
of the donor-bound electron. As is well known, a uniform 
spin polarization of the nuclei can be established through 
spin diffusion based on flip-flop processes between the 
nuclei at different distances relative to the donor cen¬ 
ter [13mi. Spin diffusion is allowed, if the energy is con¬ 
served (or nearly conserved) thereby, so that the energy 
difference between spin flips of two nuclei does not exceed 
h'yBi,, and therefore, can be compensated by the dipole- 
dipole nuclear reservoir. In our case, however, the nu¬ 
clei are exposed to an inhomogeneous Knight held Bk = 
be{r)S = feeBexp (—2r/ai), which is different for neigh¬ 
boring nuclei and is given by the electron wave function 
at the donor: 'l'^(r) = (7raf)“^ exp(—2r/a;). Here r is 
the distance from the donor center. Neighboring nuclei of 
^’’Se are separated by a distance of about R = ao/xg{^ = 
0.566/0.0758^/^ = 1.34nm [s^. The difference of the 
Knight Held at neighboring isotopes is than in the or¬ 
der of feeS'exp (—2r/aj)[l — exp (—2i?/a;)]. This should 
be compared with the local nuclear field Bl ~ 0.006 mT. 
Therefore the radial diffusion of the nuclear spin becomes 
only possible if feeS'exp(—2r/ai)[l — exp(—2B/a/)] < Bl. 
To estimate the r at which the nuclear spin diffusion be- 
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comes possible, we need to know the values for and 
S'. 

The activation energy (or donor binding energy) of the 
electron bound to the fluorine donor = 27 meV) al¬ 
lows us to estimate the localizati on rad ius of the elec¬ 
tron ai = 3.4 nm using ai = h/ with m®® = 

0.145me [30| . where me is the free electron mass in 
vacuum. This small a/ value should lead to a signifi¬ 
cant Knight field at the donor center. Using Eq. (HD 
given at the beginning of the Sec. IIV AI we obtain be = 
—370mT [d^. Then, taking into account that S = 0.07 
(see below in this section) this leads to the result, that 
the nuclear spin diffusion should be hindered within a ra¬ 
dius of about 3.9a;, resulting in this range in a spatially 
inhomogeneous nuclear spin polarization, /n,x(?')- 

Due to the spherical symmetry of the spin diffusion the 
polarization is also isotropic. Thus, the spatial distribu¬ 
tion of the nuclear spin polarization is given by: 

T (\- beS^{B - Bf^MR)exp{-2r/ai) 

“ (B - Bnmr)^ + 62^2 exp (-4r/ai) ’ ^ ^ 

We neglect here the small Bl fields. The polarized nuclei, 
in turn, act on the electrons via the Overhauser field: 

Bf^^x = ( 5 ) 

AiB5e J 

We use Eq. m to fit the experimental data for the 
induced nuclear field given in Eig. [5] The best fit is 
shown by the black curve. It has been achieved for 
be = —370mT and an average electron spin polarization 
of S' = 0.07, and reproduces all features of the experi¬ 
mental dependence. Namely, the fast changeover of the 
field sign at the resonance Snmr and the broad tails due 
to the widely extended decay. 

Eigure [10] shows schematically the behavior of the 
Knight field {Bk = foeSexp (—2r/a;)) and the nuclear 
spin polarization (/n,x) using Eq. (|4|) as a function of 
the distance from the donor. As soon as the radius be¬ 
comes larger than 3.9 a; spin diffusion becomes possible 
and equalizes the nuclear spin polarization spatially, as 
shown schematically by the red solid line showing the 
region of flat nuclear spin polarization above 3.9 a;. 

Using Eigs. 1^ and ITOl one can draw the following qual¬ 
itative conclusions: (a) close to the resonance {B = 
Snmr) the electron is exposed to a very weak nuclear 
spin polarization produced by the Knight field far from 
the donor, at the edges of the nuclear spin diffusion area. 
As the concentration of fluorine donors in this sample is 
n = 1 X 10^®cm“^, the average distance between them 
d = ( 4 ^;^) = 62 nm. This corresponds to about 

18a;, so that the donors are not located in the inhomoge¬ 
neous nuclei polarisation volumes of the neighbors, (b) 
At the center of the donor, the electrons are exposed to 
the strongest nuclear spin polarization produced by the 
maximal Knight field. This is the position far from the 
resonance on the magnetic field axis in the extended tails 
of our inhomogeneous dispersive curve in Eig.jSj A simple 



FIG. 10. (Color online) Illustration of different fields as func¬ 
tion of the distance from the donor center in units of the lo¬ 
calization radius a;. Bk, black line, is the Knight field. Here 
the simulation is done using b^S — 25.9 mT. Grey shaded area 
shows the region where the spin diffusion in the nuclear spin 
system becomes possible. In,x, red dashed line, shows the nu¬ 
clear spin polarization using B — Bnmr = b^S. Red solid line 
represents schematically the effect of nuclear spin diffusion on 
In,x, which equalizes the nuclear spin polarization spatially 
for r > 3.9 a;. Blue line shows the strength of the local fields 
Bl = 0.006 mT. 

estimation using the integral f fl/^(r)47rr^dr shows, that 
the electron is 98% confined within the volume given by 
the radius 3.9 a;. The nuclear spin polarization decreases 
at that distance from the donor by a factor of about 1200. 
If the spin diffusion border would be much closer to the 
donor or, with other words, if Bp would be bigger, the ex¬ 
tended tails would decay much faster with magnetic field. 
For example, in GaAs Bp = 0.3 mT [s^- If taking into 
account be ~ 20 mT as well as ai = 10 nm and us¬ 
ing similar considerations with S = 0.07 as for ZnSe, one 
can estimate the diffusion border to be at r « 0.3a;. This 
small borderline would give a quite homogenous nuclear 
spin polarization around the donor and lead to a nar¬ 
row NMR resonance, experimental examples are given in 
Ref. 0. 


B. Average electron spin polarization in external 
magnetic field 

The next point to clarify is the value of the average 
electron spin polarization in the range of magnetic fields 
around the NMR field. These fields are much higher 
than the half width of the Hanle curve (Bi/ 2 ). In our 
simulations using Eq. m we have estimated the aver¬ 
age electron spin value to be constant at S' = 0.07. 
The value of the average nuclear spin polarization given 
in Eq. (ED is proportional to (S • Bk) = beS^{t = 
0 )B 2 ^ 2 /(^i /2 + expected to be reduced by 

a factor of {BlBij 2 )^ ^ 1 for higher magnetic fields [l|. 
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FIG. 11. (Color online) (a) Evolution of the average spin 
components in transverse magnetic field B = {Bx,0, 0). Inset 
shows a close-np for magnetic fields from -40 to -20 mT to 
demonstrate the relative weights of the spin components, (b) 
Phase relation between the spin components. The amplitude 
of the oscillations in the red curve can be related to the red 
curve in the inset. This amplitude depends on the spin de¬ 
phasing in relation to the laser repetition period, over which 
the averaging is done. The longer the spin dephasing takes, 
the smaller the amplitude that is expected for the oscillations, 
see black curve in the inset and panel (b). 


Here {B(Bii 2 )^ ~ (20 [mT]/l [mT])^ = 400 and the aver¬ 
age electron spin transverse to the magnetic field should 
be very small due to the Larmor precession of the elec¬ 
tron spin about this field. This statement agrees with 
our experimental observations presented in Fig. [SKb) in 

Sec, mm 

However, the experimental observation of the broad 
resonance allows us to assume that there is an average 
electron spin polarization present in a wide range of mag¬ 
netic fields, which leads to polarization of the nuclei, in 
particular along the external magnetic field axis, seen as 
/N,a; in Fig. HKb). Using Eq. (13) of Ref. [2^ we can es¬ 
timate the average electron spin polarization along the y 
and z directions, where z is the optical excitation axis. 
The averaging is done over the period of the laser rep¬ 
etition, Tr. Figure fTlT a) demonstrates the evolution of 
the average electron spin components with the transverse 
magnetic field, B^ ■ The Sz component decays completely 
within a range of 10 mT (which fits well to our observa¬ 
tions in Fig. [5]), while the Sy component is slowly decay¬ 
ing up to tens of mT, see the inset to panel (a) magni fyin g 
the difference at higher fields. As mentioned in Ref. 
the vector sum of the z and y spin components decays 
with increasing transverse field B as \/B. 

Figure fTTT b) shows the corresponding evolution of the 


phase (j) between the Sz and Sy spin components. It 
allows us to conclude, that the average spin polariza¬ 
tion along the y direction should be present and have 
a constant phase shift of 90° relative to the Sz compo¬ 
nent in the range of magnetic fields above 5 mT. The 
inset demonstrates the power dependence of the phase- 
oscillation amplitude present in the simulated signal for 
different electron spin dephasing times, , in relation 
to the laser repetition period Tr. The red curve repre¬ 
sents the simulation using the realistic parameters of our 
system, T 2 = 10 ns and Tr, = 13.2 ns. 0 is the optical 
pulse area, 0 = / 2{d)E{t)dt/h, where (d) is the dipole 
transition matrix element and E{t) is the electric field of 
the laser pulse [ 2 ^. 0 = tt corresponds to the power for 
100% exciton generation. 

Therefore, the only spin component which can poten¬ 
tially provide the link to the nuclear spin polarization at 
high magnetic fields is the Sy average spin component. 
However, even if the average electron spin decays as 1/H 
the nuclear spin polarization is proportional to S^ (as 
seen shown in Sec. IHIBI) and should therefore decay as 
XjBS. This should suppress the nuclear spin polarization 
drastically with increasing magnetic field. 

On the other hand the experiments with RF field de¬ 
scribed in Sec. IHIB 31 demonstrate, that the Sy spin po¬ 
larization does not decay within the measured magnetic 
field range R = 20 — 40 mT. Also, as one can see in 
Fig. [3Kb), the amplitudes of the induced shift do not de¬ 
pend on the NMR position in the magnetic fields up to 
140 mT [ 4 ^. These two observations lead us to the sur¬ 
prising result, that the average electron spin polarization 
Sy does not change in the measured field range. 

The presented classical model gives us the possibility 
to describe the overall behavior of the signal with all its 
peculiarities. However, it requires presence of a relatively 
strong electron spin polarization Sy that does not change 
with magnetic field. Such a polarization could be caused 
by certain anisotropies intrinsic to ZnSe:F. First, the flu¬ 
orine atom itself being placed at the position of selenium 
could lead to such an anisotropy. Second, another pos¬ 
sibility could be an anisotropy of the electron spin gen¬ 
eration provided by strain in the crystal lattice. Both 
these assumptions require further investigations and will 
be presented elsewhere. 

Further possibilities to generate DNP close to RSA 
peaks may be also considered. For example, for suf¬ 
ficiently high pump power, an effective magnetic field 
along the x axis could be induced by the interaction of the 
absorption resonance with circularly polarized light (op¬ 
tical Stark effect). In that case pulsed excitation, where 
laser pulses hit the sample in phase with the electron 
spin S{t) at multiple frequencies of the laser repetition 
period Tr, induces a Sx electron spin polarization along 
the external magnetic field 4^. However, the sign of 
the induced spin polarization should depend on the rel¬ 
ative energy between the absorbtion and the excitation 
energy. This possibility was excluded by measuring the 
induced nuclear field as a function of the optical excita- 
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tion energy. The induced shift had no sign changes for 
excitation below and above the D°X resonance. 


V. CONCLUSIONS 

We have considered here that the shift of the electron 
Larmor frequencies or peaks in RSA signal in ZnSe:F 
under high excitation power is result of a dynamic nuclear 
spin polarization. This statement is confirmed by the 
position of the NMR resonance and its dependence on 
the modulation frequency /m- 

The shape and the width of the resonairce allow us to 
conclude that it should be driven by the iirhomogeneous 
Knight field acting on the nuclear system. This Knight 
field has a weak dependence on external magnetic field 
and is poiirting along the y direction. This assumption is 
confirmed by measurements, where an additional RF field 
is used to compensate the effect of the Knight field. The 
estimated values of the Knight field lead to the conclu¬ 
sion that the nuclear spin diffusion is hindered within a 


radius of about 3.9 ai froirr the fluorine donor center lead¬ 
ing in this range to inhomogeneous nuclear fields, which 
in turn contribute to the broad NMR resonance seen in 
experiment. 
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